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Abstract 

The expression for the first recoil correction to the Dirac-Coulomb spectrum is 
obtained employing the gauge invariance. 



Relativistic two-body problem in quantum electrodynamics has been exactly solved 
in the only limiting case m/M— >0, a— >0 at fixed Za (here M and m are masses of 
the constituents, Z\e\ and e are their electric charges, a = e 2 /Tic is the fine structure 
constant, H — c — 1). In this limit, an infinitely heavy nucleus holds still being the 
source of the constant in time Coulomb field. A wavefunction of the system reduces to 
that of the light particle, the electron, and obeys the Dirac equation in the Coulomb 
field, 

(dp + Pm + V c - E) if) = 0. (1) 

Expressions for the first (linear in m/M) recoil corrections to energies of the Dirac- 
Coulomb bound states were obtained several years ago by V.M. Shabaev |], 0. He used 
the perturbation theory in Za, summing up contributions of a given order in Za, linear 
in m/M. The present note is devoted to a simple derivation of the Shabaev's result, 
with only minor reference to the perturbation theory. 

As a guiding principle we will use the gauge invariance of QED. To begin with, let us 
generalize the equation ([]]) to an arbitrary gauge. Since Vc = ZaD (£\, and an infinitely 
heavy particle at rest can emit (or absorb) only zero component of the vector potential, 
we have: 

{a M (Pfj. - ZaD^o) + /3m} ip = 0, (2) 

where a = 1 by definition, p = E. The Dirac-Coulomb spectrum, that is the mutual 
arrangement of the Green's function singularities at the complex E plane[|, is certainly 
gauge-invariant . 

Now let us take into account the motion and interaction of the nucleus to first order 
in 1/M. They are described by the term 

(P- Z\e\A] 2 

(3) 

in the Hamiltonian of the system^. Here P is the operator of a nucleus momentum, 
while A is the vector potential operator acting at the nucleus site. Due to M in the 
denominator, A can be taken to act at the origin. 

In order to find the first recoil correction to an energy of the electron we need to av- 
erage the above expression over the corresponding Dirac-Coulomb eigenfunction. There 
is no problem with the vector potential operator — it emits (absorbs) photons which are 
absorbed (emitted) by the electron. Difficulties emerge when one tries to determine how 
the operator P acts on the electron wavefunction. In fact, the simple relation P = —p 
holds in the nonrelativistic limit only, when the problem is truly two-body. The rela- 
tivistic electron can propagate in both time directions so that at a fixed time slice one 
has a number of electrons and positrons with the total momentum equal to —P. 

1 D fl „'s make up the photon propagator. 

2 We consider only gauges with D^q constant in time, so that E is the integral of motion. 
3 The interaction of the electron with a nucleus proper magnetic moment, formally of the order 1/M, 
is taken into account straightforwardly, so we do not discuss it here. 
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To get rid of the problem we will start with the operator 

(ZeA) 2 

A L (4) 

whose expectation value can be easily expressed in terms of the known solution to the 
Dirac-Coulomb problem. Being only the part of the nucleus Hamiltonian (|3|), the oper- 
ator (H) is by no means gauge-invariant. This is also true for its expectation value. The 
basic idea is to reconstruct a gauge-invariant expression for the total energy correction 
from its known noninvariant part. 

Taking the expectation value of (£|) over fluctuations of the electromagnetic field we 
are left with 

&E I2 = J ^ (a,D,j{q Q )G{E - q )Dj„(-q )a„) . (5) 

Diagrammatically the right-hand side of this equation is shown in Fig.l. 



E-q 




Fig.l. "Seagull" contribution to the recoil. 



The solid line depicts G, the Green's function for the Dirac equation in the Coulomb 
field, wiggly lines represent photon propagators D^j and D Jv . As far as a gauge is not 
fixed, it is convenient to make difference between Lorentz indices corresponding to a 
nucleus interaction vertex (upper case) and those corresponding to an electron vertex 
(lower case). Overall factor 1/M allows us to take the limit M— >oo everywhere else. 
In particular, the remaining photon propagators connecting the electron line with the 
nucleus one and not shown explicitly in Fig.l, have the upper case index equal to zero 
(see (HD). As usual, Latin indices run from 1 to 3, Greek ones run from to 3. For the 
sake of brevity and later convenience, only the integral over energy flowing along the 
loop is written down explicitly. E in fl5|) is the energy of a Dirac-Coulomb eigenstate we 
average over. 

Turn now to the gauge transformation properties of (g[). The above discussion of 
the Dirac equation in the Coulomb field shows that AE^ 2 is invariant with respect to a 
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gauge transformation at the electron site, 



5D^ A = V M y? A , (6) 

where tp& are arbitrary (linear in time) functions. Hence, it remains to recover the 
invariance with respect to the transformation at the opposite 'end' of D, attached to 
the nucleus line, namely 

SD^j = Vj^. (7) 

Trying to do this we cannot use components of D with a spatial upper case index. 
Actually, those components of D emerge in an expression for the energy correction due 
to the operator A acting at the nucleus site. Since the quadratic in A effect is already 
taken into account by (fD, the only thing that can help us reads 

« M -D M o = Aio - ««Ao- (8) 

It does not spoil the invariance with respect to (^). On the other hand, its gauge 
variation with respect to the evident extension of (0), 

SD^ (qo) = iq (f^ (9) 

does compensate that of a^D^j in the linear combination 

a^^V^o) + — <V l^J> D M<lo)} ■ ( 10 ) 



By this means the gauge-invariant expression for the recoil correction takes the form: 

AE = ^[^ (11) 
M J 2tt k j 

(a^ (d^o) + hi [Vj, i^o(go)A G(E - g ) ^D Ju {-q ) - h [Vj, D 0u (-q )]\ a v \ . 

Recall that «o = 1. Unfortunately, this expression is meaningless until we define how to 
treat the new singularity at q = 0. 

As soon as the gauge invariance is maintained we can choose the mostly convenient 
gauge. Without question this is the Coulomb one. Then the total energy shift ( pT|) is 
naturally broken up into four terms: 

AE = AE CC + &Ecm + A£ MC + AE MM , (12) 

the last of which is nothing but the Coulomb gauge version of A£^ 2 (see (||) and Fig.l), 
i. e. the double magnetic exchange contribution to the energy shift. The third and the 
second terms comprise the correction arising due to a single magnetic exchange. Their 
origin at (0) is the term — Z\e\(PA + AP)/2M . Finally, the pure Coulomb contribution 
AEcc is just the mean value of the nucleus 'kinetic energy' P 2 /2M. 
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To find a prescription according to which the 1/go-singularity should be passed by 
one can exploit its independence of Za and analyze the corresponding expression pertur- 
batively in Za, to the lowest nontrivial order. Very natural result for the single-magnetic 
exchange reads 

1 1/1 1 

g 2 \q -iO g + iO 

For example, the third term in ( |12"D can be represented diagrammatically by the sum of 
two graphs shown in Fig. 2. There the thick line depicts the propagator (±go + ?0) _1 of 
the infinitely heavy nucleus having the energy M ± g , while the dashed line shows the 
interaction through the Coulomb electric field Z\e\ [V, D 00 \. 




E-q E-q 




M + q M - g 

Fig. 2. Single magnetic exchange. 



Simple perturbative analysis of the pure Coulomb contribution AE C c suggests that 
corresponding prescription has the form: 

11/1 1 
I^2l v (go- i O) 2 + (go + i O) 2 

Actually, here we have the sum of the nucleus propagator derivatives resulting from the 
expansion of (P 2 /2M ± go + iO)^ 1 . To obtain Feynman diagrams for the pure Coulomb 
contribution, one need only substitute wiggly lines in Fig. 2 for dashed ones as well as 
the nucleus propagator for its square. 

Now, when the integral in (|TTD is completely defined one can easily check that this 
expression equals the sum of recoil corrections found by Shabaev in the more straight- 
forward way [[l], (in 0, the overall sign of the expression for AE MM is corrected). 

Recall that our starting point was the expectation value (0). Its perturbative ex- 
pansion can be readily appreciated to begin with (Za) 5 . To be certain that the corre- 
sponding expansion for the total correction (|ll|) begins with (Za) 2 , let us consider in 
greater detail the pure Coulomb contribution which alone survives the transition to the 
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nonrelativistic limit. As a byproduct we will see how the expectation value \P 2 /2Mj 
looks in terms of the solution to the Dirac-Coulomb problem. 

Evaluating the integral with respect to q in AE C c according to the prescription (|14D 
together with the standard rules for the Dirac-Coulomb Green's function one readily 
obtains 

&Ecc=^{p(A+-A-)p), (15) 

where A± are the projection operators to sets of positive- and negative-energy Dirac- 
Coulomb eigenstates correspondingly. Passing from (|TT| ) to (|15D we used the Dirac 



equation ([I]). In the nonrelativistic limit A + — >1, A_^0 and ( JIB] ) reduces to the well- 
known result, 
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